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Abstract 

The objective of this study was to identify barley leaf proteins differentially regulated in response to drought and heat 
and the combined stresses in context of the morphological and physiological changes that also occur. The Syrian 
landrace Arta and the Australian cultivar Keel were subjected to drought, high temperature, or a combination of both 
treatments starting at heading. Changes in the leaf proteome were identified using differential gel electrophoresis 
and mass spectrometry. The drought treatment caused strong reductions of biomass and yield, while photosynthetic 
performance and the proteome were not significantly changed. In contrast, the heat treatment and the combination 
of heat and drought reduced photosynthetic performance and caused changes of the leaf proteome. The proteomic 
analysis identified 99 protein spots differentially regulated in response to heat treatment, 14 of which were regulated 
in a genotype-specific manner. Differentially regulated proteins predominantly had functions in photosynthesis, but 
also in detoxification, energy metabolism, and protein biosynthesis. The analysis indicated that de novo protein bio- 
synthesis, protein quality control mediated by chaperones and proteases, and the use of alternative energy resources, 
i.e. glycolysis, play important roles in adaptation to heat stress. In addition, genetic variation identified in the pro- 
teome, in plant growth and photosynthetic performance in response to drought and heat represent stress adaption 
mechanisms to be exploited in future crop breeding efforts. 

Keywords: Abiotic stress, barley, drought, heat, proteomics, Rubisco activase, yield. 



Introduction 

Drought and heat are among the main abiotic stresses dra- 
matically limiting crop growth and productivity worldwide 
(Wang et al., 2003). In the field, the co-occurrence of several 
abiotic stresses, rather than an individual stress condition is 
most damaging to crop production (Mittler, 2006). For exam- 
ple, the combined effects of heat and drought on yield are 
more detrimental than the effects of each stress alone, as seen 
in sorghum (Sorghum bicolor L.; Craufurd et al, 2008), wheat 
(Prasad et al, 2011), and barley (Savin and Nicolas, 1996). 
However, most studies to date have addressed the effects of 



single stresses on plant performance (Ugarte et ah, 2007, 
Harb et ah, 2010), and little is known about the molecular 
mechanisms underlying the acclimation of plants to a com- 
bination of different stresses (Mittler, 2006). Recent studies 
have revealed that the response of plants to a combination 
of different abiotic stresses is unique and cannot be directly 
extrapolated from the response of plants to each of the differ- 
ent stresses applied individually (Rizhsky et al, 2004; Ahmed 
et al., 2013). Breeding of stress-tolerant crops is the most effi- 
cient strategy to maintain yield in stress-prone marginal land. 



© The Author [201 3]. Published by Oxford University Press [on behalf of the Society for Experimental Biology]. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/ 
by-nc/3.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, 
please contact journals.permissions@oup.com 



3202 Rollins et al. 



It is thus important to identify genetic resources with high 
stress tolerance and to understand the mechanisms contribut- 
ing to adaptation to stresses typically co-occurring in the field 
such as heat and drought. 

Barley (Hordeum vulgcire L.) is worldwide among the most 
widely cultivated crops in marginal environments and is often 
the most common crop in the driest rain-fed farming areas as 
it is well adapted to abiotic stresses (Baum et al, 2007). Barley 
has thus been selected or bred for specific adaptation to abi- 
otic stresses in geographically distinct areas of the world. 
This adaptation of genetically diverse germplasm to similar 
environmental conditions over a wide geographical range 
can be exploited for breeding and germplasm exchange. For 
example, barley germplasm bred by the International Center 
for Agricultural Research in the Dry Areas, Syria (ICARDA) 
for the marginal environments of West Asia and Northern 
Africa (WANA) showed good adaptation to the dry Southern 
Australian environments (Coventry et al, 2004) and vice versa 
(Shakhatreh et al, 2001). This germplasm is also interesting 
in order to study the genetic basis of adaptation to stress 
in genetically diverse genotypes. Stress adaptation in barley 
has been attributed to the plasticity of morphological traits 
such as biomass production, plant growth, tiller number, and 
peduncle extrusion (von Korff et al, 2008; Shakhatreh et al, 
2010). Growth-related responses to stress are typical of stress 
avoidance strategies which allow the plant to maintain a 
homeostasis despite changes in the environments. Barley has 
also been characterized for genetic variation in physiological 
parameters such as relative water content and chlorophyll flu- 
orescence parameters under stress (Oukarroum et al, 2007; 
Ahmed et al, 2013). These differences are indicative of differ- 
ences in stress-tolerance mechanisms which allow the plant to 
maintain cellular activities under stress (Bartels and Sunkar, 
2005). Studies on molecular changes in response to stress 
in barley have primarily relied on quantification of mRNA 
changes under stress (Talame et al, 2007; Guo et al, 2009 
von Korff et al, 2009). These studies showed that a large pro- 
portion of the transcriptome responded to drought includ- 
ing genes implicated in stress signalling and stress response. 
However, while mRNA transcript levels can dictate protein 
abundance, the differential expression of the two macromol- 
ecules is not always well correlated (Stylianou et al, 2008; 
Baerenfaller et al, 2012). Proteomic approaches provide 
information missing in DNA or mRNA analysis methods 
in that they focus on the actively translated portion of the 
genome. Stress resistance is conferred by the proteins, which 
function in stress signalling, transcription regulation, cellular 
detoxification, protection of macromolecules, and a panoply 
of other processes. In recent years, methods for the analysis 
of the proteome have advanced considerably, and together 
with emerging sequence information in model crops like 
barley (The International Barley Sequencing Consortium, 
2012), the plant proteome has become increasingly useful 
for understanding gene function and networks in response 
to environmental stimuli. Conventional gel-based proteom- 
ics has proved useful in barley research to quantify changes 
in protein abundance in grains during development (Finnie 
et al, 2006), in roots in response to salt stress (Witzel et al, 



2009) and in shoots in response to heat stress (Siile et al, 
2004). Proteomic research using fluorescent labels and two- 
dimensional difference gel electrophoresis (2D-DIGE) has 
been successfully applied in barley for the identification of 
proteins associated with malting quality (March et al, 2012) 
and in Arabidopsis (Shi et al, 2011) and wheat (Gao et al, 
2011) to identify proteins responsive to salt. By quantifying 
changes in protein abundance, one can gain insight into the 
biochemical processes that underlie the plant's morphological 
and physiological acclimations to abiotic stress. 

The objectives of this study were to study the molecular 
basis of abiotic stress responses in two drought-adapted, 
genetically diverse genotypes, the Syrian landrace selection 
Arta and the Australian barley cultivar Keel. The specific 
aims were to: (i) characterize the physiological and mor- 
phological responses to drought and heat stress; and (ii) 
identify leaf proteins differentially regulated in response to 
drought and heat stress using a proteomics approach based 
on DIGE and mass spectrometry. The changes in protein 
abundance are discussed in the context of the physiologi- 
cal and morphological trait plasticity that occurs due to the 
abiotic stresses. 

Materials and methods 

Plant material 

Two cultivated barley two-row genotypes, Arta and Keel, were 
grown in a growth chamber to evaluate their performance. Arta is a 
pure line selection from a Syrian landrace adapted to the driest areas 
in Syria and is winter barley. These genotypes were selected because 
they are genetically diverse but adapted to similar drought-prone 
environments, thus allowing the study of the diversity of molecular 
and phenotypic changes in response to drought and heat. 

Drought and heat treatments 

Response to the single or combined effects of drought and heat 
treatments applied at the generative stage was tested in a duplicated 
growth chamber experiment. Plants were sown in 96-well trays, strat- 
ified at 4 °C for 4 d, and then grown in 8/16 light/dark short-day 
(SD) conditions for 24 d. Plants were transferred to 16/8 light/dark 
long-day (LD) conditions for 2 d to acclimate before being potted in 
4 1 pots containing 1.8 kg soil with three plants in each pot. Plants 
remained in LD for the remainder of the experiment. The SD fol- 
lowed by LD treatment was applied to synchronize flowering. The 
relative humidity of the chambers was set to 50%, the light inten- 
sity was 350 umol photons mT 2 s _1 and the temperature was set to 
21/17 °C light/dark. The field capacity (FC) was calculated as the dif- 
ference in weight between fully hydrated soil and dried soil (Colman, 
1947). The soil water content (SWC) of potted plants was adjusted to 
50% of the FC. For the drought treatment, the SWC was reduced to 
1 5% FC by controlled withholding of water; pots were weighed daily 
and watered to match the weight of the heaviest pot. The reduction 
of the SWC was equal across all drought-treated pots. The SWC of 
1 5% FC for the stressed plants and 50% FC for the control plants was 
maintained until plant maturity. To control for the added weight of 
the growing plants, the volumetric water content of random pots was 
checked weekly using a TDR 100 soil moisture meter (FieldScout) 
fitted with 12 cm probe rods. When the SWC reached 15% FC in the 
drought-treated plants, the heat treatment was applied to a subset 
of the well-watered and drought-treated plants by moving the pots 
at ZT 3 to an identical chamber set to 21 °C and gradually raising 
the temperature to 36 °C over 4h. Heat-treated plants remained at 
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36/32 °C light/dark for 1 week, at which point the temperature was 
decreased to 21 °C over 4h. Detailed information on phenotyping 
can be found in Table 1 and in the Supplementary Materials and 
Methods (available at JXB online). 

Two-dimensional gel electrophoresis 

Leaf proteins from three biological replicates were extracted using a 
trichloric acetic acid/acetone precipitation modified from Cascardo 
et al. (2001) and fluorescently labelled as detailed in Supplementary 
Table SI and Supplementary Materials and Methods. The protein 
sample was diluted to a total volume of 340 ul in rehydration buffer 
consisting of 7M urea, 2M thiourea, 2% CHAPS, 20 mM DTT, and 
0.5% biolyte 3-10 ampholytes (Bio-Rad) and applied to 18 cm immo- 
biline strips pH 3-10 NL (GE Healthcare Life Sciences). Strips for 
DIGE gels contained a total of 1 50 |ig protein while strips for gels 
to be post-stained with PageBlue (Fermentas life sciences) contained 
500 (.Lg protein. Focusing was accomplished at 20 °C using a Protean 
isoelectric focusing cell (Bio-Rad) with the following conditions: 
14h passive rehydration, 250V for 15min, ramping to 2000V for 1 h 
45min, and ramping to lOkV for 3h, before maintaining the volt- 
age at lOkV for a total of 30kV/h. The resulting strips were equili- 
brated in 2D equilibrium buffer (0.1 M Tris, 6M urea, 30% glycerol, 
2% SDS) containing 2% DTT for 15min and then with 2D equi- 
librium buffer containing 2.5% iodoacetamide (Sigma) for 15min. 
Equilibrated strips were placed on 1 mm thick 12% SDS-PAGE gels 
sized 26 x 20 cm and covered with 0.5% agarose (Bio-Budget). The 
second dimension was separated using 12 mA/gel for 12 h with the 
EttanDaltSix electrophoresis system (GE Healthcare Life Sciences). 

Protein identification 

Protein gels were imaged and spots were detected as detailed in 
Supplementary Materials and Methods. Spots from 2D gels were 



Table 1 . Summary of traits measured in the genotypes Arta 
and Keel grown under control, drought, heat, and combination 
treatments 



Trait 


Abbreviation 


Unit 


Grain yield 


GY 


g 


Total biomass 


BM 


g 


Harvest index 


HI 


gg~ 1 


Plant height 


PH 


cm 


Peduncle extrusion 


Pedex 


cm 


Spike number 


SN 




No. of aborted spikes 


AS 




Grains per spike 


GS 




Thousand kernel weight 


TKW 


g 


Days to maturity 


DM 


days 


Total water used per pot 


WU 


I 


Water use efficiency of grain yield 


WUE 


gr 1 


Leaf temperature on day 1 


LT_1 


°c 


Leaf temperature on day 3 


LT_3 


: c 


Leaf temperature on day 7 


LT_7 


°c 


Relative water content on day 1 


RWC_1 


% 


Relative water content on day 3 


RWC_3 


% 


Relative water content on day 7 


RWC_7 


% 


Maximum PSII quantum yieldat day 1 


Fv/Fm_1 


Arbitrary 


Maximum PSII quantum yield at day 3 


Fv/Fm_3 


Arbitrary 


Maximum PSII quantum yield at day 7 


Fv/Fm_7 


Arbitrary 


PSII performance index at day 1 


Pl_1 


Arbitrary 


PSII performance index at day 3 


Pl_3 


Arbitrary 


PSII performance index at day 7 


Pl_7 


Arbitrary 



excised using the Proteineer sp II and tryptically digested and spotted 
on MTP 600/384 AnchorChip plates using the Proteineer dp system 
(Bruker Daltonics). Aliquots of the digests were automatically spotted 
on MTP 600/384 AnchorChip plates for subsequent mass spectromet- 
ric analysis according to the thin-layer protocol of Gobom et ah (2001 ) 
using the Proteineer dp robot. Peptide mass fingerprints were obtained 
using the Ultraflex III MALDI ToF/ToF mass spectrometer (Bruker 
Daltonics). The resulting spectra were processed into peak files with 
the flexAnalysis version 2.4 software (Bruker Daltonics) by means of 
the sophisticated numerical annotation procedure algorithm. Peak 
data were imported into the ProteinScape database system version 3.0 
(Bruker), which initiated Mascot version 2.3 (Matrix Science) searches 
against the UniProt knowledgebase (http://www.uniprot.org/) for 
H. vulgare (release 2011_06) and the DFCI Barley Gene Index (http:// 
compbio.dfci.harvard.edu/tgi/) version 12 genome database. 

Data analysis 

Analysis of variance (ANOVA) and correlations were performed 
using SAS version 9.1 (SAS Institute) and are described in details in 
Supplementary Materials and Methods. The regulation factor, the 
log2-transformed fold change in the normalized spot intensity, was 
calculated between 36 °C and 21 °C treated plants across genotypes 
and drought treatments as well as between Arta and Keel across both 
stress treatments. Gene ontology terms for identified proteins were 
retrieved from AgBase (McCarthy et ah, 2006). Singular enrichment 
analysis was performed with the AgriGO toolkit (Du et ah, 2010) as 
outlined in Supplementary Materials and Methods. 

Results 

Comparison of effects of drought, heat, and 
combination treatments on plant performance 

Morphological changes caused by the treatments are depicted 
in the representative photographs of control and stress-treated 
Arta and Keel plants taken the end of the heat treatment, 
9 days after heading initiated (Fig. 1). Drought caused strong 
reductions in plant height as compared to control in Arta, but 
not in Keel. In addition, the drought, heat, and combination 
treatments caused noticeably stronger senescence of lower 
leaves in Arta than in Keel. 

Overall, the drought treatment had a stronger effect on 
morphological traits than on physiological traits. For exam- 
ple, the water regime explained 26% of the variation in 
grain yield (GY), 57%> of the variation for biomass (BM), 
and 79%o of the variation for spike number (SN), while it 
explained only up to 18% of the variation in relative water 
content (RWC) and less than 1%> of variation in photosyn- 
thetic performance index (PI) (Supplementary Table S2). 
In contrast, the heat treatment had significantly stronger 
effects on physiological traits than the drought treatment; 
it explained up to 54%o of the variation in water use effi- 
ciency (WUE), 34%o of the variation in RWC, and 74%o of 
the variation in PI. The heat treatment did not have strong 
effects on plant growth; it only explained 8%o of the vari- 
ation in BM and 2%o for plant height (PH). The strongest 
phenotypic changes were observed under the combination 
treatment (Table 2). The combination treatment caused the 
strongest reduction in GY. GY under the combination treat- 
ment was significantly lower in Keel than GY under heat or 
drought alone, while significant reductions in GY caused by 
drought and heat were not significantly different from each 
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other (Table 2). While the reductions in GY under drought 
and heat were similar, drought and heat had different effects 
on yield component traits. Whereas the drought and combi- 
nation treatments caused significant reductions in BM and 
PH, heat did not significantly affect these traits. In addition, 
the spike number (SN) was significantly reduced by drought 
and combination treatments as compared to controls, but 
not significantly affected by heat in both genotypes. In 
contrast, heat and combination treatments significantly 
increased the number of aborted spikes (AS) and decreased 
thousand kernel weight (TKW) in Arta and Keel, while 
drought did not have significant effects on AS and TKW as 
compared to control conditions. The harvest index (HI) was 
significantly reduced under drought, heat, and combination 
treatments as compared to controls, while the reductions in 
HI were also significantly greater under heat and combina- 
tion treatments as compared to the drought treatment in 
both genotypes. The number of grains per spike (GS) and 
peduncle extrusion (Pedex) were reduced under drought, 
heat, and combination treatments, but the differences were 
not significant. The heat treatment had thus strong effects 




Keel:2il°-C 




Fig. 1. Representative Arta and Keel plants after 7 days of control 
(15%) soil water content. Each 4 I pot contains three plants. 



on generative traits; it explained 32% of the variation of AS, 
and 72% of variation in TKW. In contrast, drought showed 
the strongest effects on vegetative traits such as vegetative 
biomass and plant height. 

Under drought, GY was positively correlated to BM (0.70) 
and PH (0.56), while under heat and combination treatments 
GY correlated with BM (0.57 and 0.46, respectively), but 
not with PH (Supplementary Table S3). Under heat alone, 
GY correlated negatively with AS (-0.54) and positively with 
TKW (0.58). Different yield components were thus correlated 
to GY under the different treatments. 

WUE, RWC, maximum quantum efficiency of photosys- 
tem II (Fv/Fm), PI, and leaf temperature (LT) were signifi- 
cantly affected by heat and combination treatments, but not 
by drought (Table 2). WUE dropped from 1.6 and 1.7 under 
control to 0.7 and 1 .0 under heat and 0.9 and 0.9 under com- 
bination treatments in Arta and Keel, respectively. In addi- 
tion, at 3 and 7 days after start of the treatment, Arta and 
Keel plants had significantly lower RWC due to the com- 
bination treatment (Table 2). The integrated chlorophyll 
fluorescence measurement PI was significantly reduced 




Keel 36°C 




°C) or high (36 °C) temperature and control (50%) or drought 



Leaf proteome alterations in responses to drought and heat stress in barley | 3205 











































LO 


UO 


CM 
















CO 
















CO 
















CO 


CO 


CM 










x 


CO 


■ 




■ 


■ 


■ 


■ 


■ 




0 


■ 


■ 


■ 




■ 


■ 




CO 






Jd 




■ 






CO 


^~ 


CO 






1" 










0 


CO 






















CM 














LO 


CM 






CM 








CO 


CO 


CO 




CO 


CO 














\j 

CD 
CO 










































CD 


CM 














CM 














































LO 




CM 




I s - 


fs_. 




■ 


- 




CO 
- 


- 


- 


- 


in 


- 


t — 1 
- 




- 


CO 
- 


- 


P 

- 


LO 


- 


c 

™ 




LO 


O 


OO 




LO 


0 




CO 


CO 


CM 


0 


LO 


LO 


■ 

CO 


CO 




■ 

0 


0 


0 


O 




■ 

O 


0 












°f 








































£ 










































^ 


^ 








Ui 
3 


c 


u 


13 


-0 
O 


§ 


a 


■O 




CB 


13 






(J 




■o 


13 


S 


"8 




LO 


CM 
CO 


■st 

CO 


-0 


U 


13 






OO 


CO 


CO 


CO 


CO 


co 


CO 


O 


CO 


cn 


CM 


0? 


CM 


0 


CM 


CNj 


CO 






h~ 


LO 




O 


CO 


£ 


CD 


CM 


05 


0 


LO 


1 


■vt 


CO 


d 


CD 


O 


CO 


0 




h- 


h- 


G} 


CM 


d 


O 


0 


O 


CM 




0 


Q 








-st 


1— 






CM 


03 






CO 


CO 


CO 






LO 


























































O 
































































O 


0 




h- 


CM 


r>- 




■st 




CM 


CO 


CO 


0 


CO 


h- 


CM 


^ 


UO 


CO 


CO 


00 


CO 


UO 


CO 


CO 


CO 


CD 


CO 


d 


CO 




LO 


co 






d 




d 


o5 


d 


cj 


CO 


d 


d 


d 


d 


d 


d 










1— 




■st 


1 


CM 






CO 


1- 






CO 


CO 


CO 


CJ> 


CO 




































































O 








CO 
















CO 
















uo 














CO 


CO 


0 


I s - 


CO 


r-- 


S- 


CM 


CO 




CO 


CM 


CO 


LO 


03 




00 


, 






co 


uo 




^ 


co 


c 


d 


d 


d 




id 






d 




CO 




d 


d 


d 


CM 


LO 




d 


d 


d 


d 




d 




CO 


jg 








CO 


l 








CM 


00 






CO 


CO 


CO 


LO 


uo 


>* 














3: 




















































CJ> 


c 


■g 


"p 




■g 


■Q 


u 


13 














3d 


■g 






"g 




CO 










3 

£ 


CO 


CM 


CO 


CM 


CO 


to 


to 


CO 


CO 


C53 




CD 


d 


LO 




■st 




CO 


CO 






r-- 


0 


CO 




CD 


CO 


CM 


O 


CM 


l6 


d 


CO 


d 


d 


b 


03 


d 


cb 


id 


d 






d 


d 


d 


d 


d 






Q 










Sf 


1 








CM 


03 






CO 


CO 


CO 




CO 


CO 






















































uo 


CM 


CM 
















CO 
















LO 






























x 


0 


CO 


LO 


O 


-st 


0 


0 


0 


'st" 










CO 


CO 


<M 


cn 


1^ 


CO 


CO 


00 


CM 


00 


CO 


eel 


CO 


00 


CO 


0 






LO 


CO 


CO 


CO 


CO 


CO 


CO 


LO 


LO 




CO 


CO 




0 


d 


d 




d 


d 










LO 




CM 








35 






CM 


CM 


CM 


o> 


O) 


00 


































































O 

0 

T - 








LO 




CO 




























0 


CO 


0 










O 


CO 


'st 


h- 


I s - 


0 


CD 


CM 




O) 


CO 


LO 


I s - 


■st" 


CN 


m 


O 


CO 




00 


CD 


CM 


CO 


c 


CD 


CM 


0 


LO 




CM 


d 


O) 




CO 


CM 




d 


d 


d 


O 


uo 


CO 


0 


d 


d 


CM 


CM 


CM 


CM 


i 










T 








"sf 


CD 














CO 


LO 














+j 

.C 

O) 


c 


-Q 


-Q 


U 






■H 










_ 

CM 












abed 


Id 


cb 

CM 




-8 

LO 








3 

£ 


co 






LO 


OO 


CO 


CO 




?0 


UO 




'sf 


CO 


0 


cb 


00 


T- 


CO 




CO 


CO 


CO 




cb 


CNJ 


CO 


rsj 


LO 


d 


d 


LO 




CM 




■st 


CO 


CO 


CM 


CM 


CM 




LO 


00 


LO 


d 


d 


d 


CO 


d 


d 


Q 










■st 


| 








■st 


CO 






CM 


CM 


CM 


CO 


1^ 


r^. 






















































CM 


CO 


CO 








TO 


Max 


LO 




CM 
CO 




LO 




CO 






LO 


CO 
CM 


CO 










CO 




CO 


CO 


00 




r — 


CO 


Art 


CO 


d 


d 


CM 
'st 


CM 
1 


d 


LO 


d 


LO 


O 


■st 


CM 


LO 
CM 


d 

CM 


LO 
CM 


CO 
CO 


d 

CO 


UO 
CO 


d 


d 


d 




d 










































CM 


CO 


CO 








O 








CO 
















CO 


















0 


■ 
















CM 


















LO 




CO 






^ 


r" 










CO 


T - 


c 


CO 


■ 


d 


d 


d 


CO 


■ 


if) 


d 








d 


d 


00 


LO 


CO 


CO 




0 


0 


• 


d 


CM 


<N 


I 








CO 


1 








CO 
















CO 


uo 














£ 

CJ> 


c 


13 




u 




-Q 




y 


















■Q 


u 




(D 




0 








3 


co 


■& 


cb 


4 


d 




"d 


CM 


CM 


d 




■st 


00 


CO 


'd 


CO 




So 


0- 


CO 




00 


d 


d 


CD 


£ 


<D 


CD 


CO 


d 


03 


9 




d 




d 


d 


03 


T — " 


CM 


d 


CM 




d 


d 


d 


d 


d 


CO 


d 


d 


Q 










CO 


7 










03 






CM 


CM 


CM 


CO 


























































CM 


LO 


CM 
















LO 
















CD 
















o> 


CO 


00 










X 


CM 


CO 


'st 


CO 


00 


CO 


0 




0 




O 


CO 


CO 


03 


LO 


O 


cq 


CM 


t-- 




I s - 




CO 


00 


"5 


co 


CO 


CO 


O 


d 




d 


d 




d 




CO 


1— 


LO 




1^ 


0 




d 


O 


d 


d 


CM 


CM 


1 


a> 






CM 




LO 




CO 






CO 


0 






CO 


CO 


CO 


0 


00 


CO 






















































1— 




CD 








p 








h- 




































CO 








c 


CO 


CO 


CM 


*st 


O 




CO 


CM 


0 




CO 


CO 


CO 


03 


CM 




10 


CO 


r — 




LO 


CO 


CM 


CO 


CD 


i 


































d 
















CO 




1_ 




"sT 


' 


CM 






CM 


CSJ 






CO 


CTJ 


VJ 


CiJ 




cu 














:rol 


c 




-Q 


T3 
■O 




-Q 




■0 




■g 




So 








a 




abc 


-0 


i3 
5; 


u 
0 


to 
00 








c 


co 


% 


CO 


CO 


00 


OO 








■sf 




CD 


p 


CM 


03 


CD 


CO 


LO 








I s - 


CO 


00 




0 


0 


CD 


co 


d 


d 


LO 


d 


cb 


d 


d 


00 


cb 




CO 


ib 


■»t 


CO 


d 




d 


d 


d 


c\i 






O 










LO 


1 


CM 






CO 


0) 






CO 


CO 


CO 


r^- 


CO 


























































CO 


03 
















CM 














































X 


CO 


CO 


LO 


CM 


03 


CO 


0 




CM 




d 




CM 


CO 


CD 


CN 






CO 


CO 


00 


CO 


CO 






co 




CO 


d 




CM 


d 


CM 


LO 




LO 


d 




d 


d 










d 


d 


d 


d 


d 


d 


rta 






CM 




LO 




CO 


CM 






O 






CO 


CO 


CO 


5 


co 


N- 




















































CO 


0 


CO 
















































CO 


CM 










p 




to 


CO 




O) 


CO 


CO 


O 


CO 


T— 




CM 




00 


r-- 


CO 


CO 


CO 


CM 






iS- 


uo 


CO 


CD 


c 


oj 


LO 


d 


d 


CO 


d 




d 


CM 


CM 


d 


d 


CM 


d 




LO 


d 


CO 


d 


d 


d 


'- 


d 


d 


CD 


i 










| 






















LO 


















CO 












































« 








:rol 


c 


"8 


-Q 


















LO 










abc 


bed 








S 








c 


co 


"co 




CM 


bo 


I s - 


d 


i- 




LO 




CO 




q 


q 


CM 


CO 


CM 


% 






N 


d 


CO 


"cb 


0 


<D 


LO 


CO 


d 




co 


LO 




d 


d 






d 


LO 


id 


d 






d 


d 


d 


d 








O 












1 


CM 






CO 


O) 






CO 


CO 


CO 


CO 


s. 
























































CO 




0 
















OO 
















CD 


















CM 












X 


CM 


T- 


h- 


<st" 


OO 


O 




CM 


O 








0 


CO 


-st 


■<t 


CO 




CO 


CO 


00 


CO 


LO 




■5 


co 


CO 


CO 


0 


cb 




CO 


C\j 




ib 


03 




CM 


LO 


'st 


CO 


CO 




d 


0 


d 


d 




d 


d 


a> 






CM 




LO 


T 


CO 






LO 


O 
1— 






CM 


CM 


CM 


O) 




CO 
























































CO 


LO 














^t 
















I— 


















0) 


CO 








p 




CM 


LO 


■st 


LO 






0 


CO 








CO 




O 


CO 


CO 






CO 




I s - 




CD 


CM 


c 


*st 


OO 


d 


cb 


d 


d 


d 




CO 


CO 


-st 


d 


d 


d 


06 


LO 


LO 


d 


d 


d 


d 


CM 


CM 


d 


T - 

CM 


I 








<=t 


1 








CO 


CD 








CM 




CD 


1^ 


CO 














£ 










































■Q 


-0 










san 










?- 








■Q 




O 












■Q 




LO 
CO 


CD 


eo 




0 




Jon 


I s - 


CM 




LO 




CM 
















d 


^t 


CNI 


LO 








CO 


&5 




00 






O 


d 


CM 


LO 

1 


CM 


d 


CM 


, st 


00 


cb 










CO 




d 


d 


d 


0 


CO 


d 


CM 






CM 




LO 


CM 






■st 


CO 




CM 


CM 


CM 


CO 


CO 


















































CO 


00 


CD 
















CO 
































CO 














X 


1 ~; 


CO 




CO 


CO 


O 




CM 


CO 




CM 




LO 




■st 


CO 


CO 


CM 


CO 


CO 


CO 


CM 






rta 


co 


LO 




d 


d 


C\j 


d 


d 




d 




d 


CM 




^t 




CO 


CO 


d 


d 


d 


d 


16 


d 


CO 






CM 




LO 


1 


CO 






CO 


0 






CM 


CM 


CM 


CO 


CO 


05 














< 








































CO 


0 


CO 








p 








0 
















C33 


















CO 










c 


LO 


CM 




0 


CD 


0 


CO 




CD 




CO 


00 






CO 


CM 






5 


r~- 


CO 




03 


CM 


1— 


I 


CO 


CO 


0 






CO 


0 


0 


O 




LO 


O 




CO 






d 




0 


0 


d 


CM 




CM 


CM 








'st 


1 










CO 












CO 


r~- 
















£ 








































CO 


0 


u^ 

CO 










c 


CO 




cb 






u 


0 










13 














CO 


CM 


CM 








c 


co 


CO 




LO 


cb 




cb 


LO 


CM 


CM 


_ 


LO 


CO 


CO 


,[ d 


d 


CO 


uo 


d 


CO 


CO 


CO 


00 


LO 


O 


0 


<D 




d 


d 


d 


CO 


CM 


d 


CM 


d 








d 


i— 


d 


d 


DO 


CM 


d 


d 


d 


CO 


d 


d 


0 






CM 




"st 


1 


CM 






LO 


03 






CM 


CM 


CM 


CO 


CO 


CO 


















































CO 






CO 

1 


N- 














































1 




















8 














LU 




CO 




^1 
0 


1 

0 


0 


E 
u_ 


E 
u_ 


E 

Li_ 




CO 

1 




'is 




>■ 




— 


X 


"O 

© 


z 


cz> 


00 






WU 


WU 






b' 


g 




g 


L? 






'"l 


1 


1= 




0 


CD 




CL 




CO 


< 


0 


H 


O 








LX 




LX 








Q_ 


D_ 


Q_ 



3206 Rollins et al. 



under heat and combination treatments in both genotypes 
at all analysed time points. In contrast, Fv/Fm was only 
significantly reduced in Keel under heat and combination 
treatments at all three time points analysed. LT (as opposed 
to the ambient temperature) was significantly increased by 
drought 1 day after treatment start in Arta plants (22.8 °C) 
compared to well-watered controls (20.3 °C), but no sig- 
nificant differences in LT were detected between drought- 
treated and well-watered Keel plants (Table 2). LT was 
significantly higher in Arta and Keel plants grown under 
heat than plants grown under control temperatures, with 
mean values ranging from 33.2 to 37.0 °C and from 20.3 to 
23.0 °C, respectively. After 7 days of heat treatment, both 
Arta and Keel plants experienced higher LT when addition- 
ally treated with drought (36.4 and 37.2 °C, respectively) as 
compared to the heat treatment alone (33.2 and 34.6 °C, 
respectively). 

The experiment revealed the strongest genetic varia- 
tion for the traits PH (24%), AS (17%), and Fv/Fm (15%) 
(Supplementary Table S2). Arta showed a significantly 
stronger reduction in PH under drought than Keel. In addi- 
tion, heat resulted in a significantly higher amount of AS in 
Arta than in Keel (Table 2). Finally, in Keel, the combination 
treatment resulted in a significantly lower Fv/Fm and PI as 
compared to the heat treatment alone and as compared to the 
Fv/Fm and PI under the combination treatment in Arta. The 



two genotypes thus showed different physiological responses 
to the heat and combination treatments. 

Characterization and quantification of the barley leaf 
proteome 

Spot detection on Coomassie-stained gels (Fig. 2) revealed a 
total of 525 distinct spots that were submitted for identification 
by peptide mass finger printing and subsequent peptide frag- 
mentation analysis. Mascot search against the Uniprot data- 
base for H. vulgare or the DFCI Barley Gene Index database 
revealed the identity of 296 proteins from a total of 525 spots. 
The 296 identified spots were matched to a total of 145 unique 
accessions. Spots identified with the same protein accession 
were considered to be isoforms. Sixty-two of the proteins iden- 
tified in the barley leaf proteome had isoforms present. 

Using a singular enrichment analysis with gene ontology 
terms and the UniProt database, over half of the identified pro- 
teins were present in intracellular compartments (59%) and the 
cytoplasm (52%), as well as in organelles (51%) and were signifi- 
cantly enriched compared to the background (22, 9, and 1 5%, 
respectively; Supplementary Fig. SI). The largest fold enrich- 
ment was for the term 'plastid', which was present in 45% of 
leaf proteins and in less than 2% of background proteins. 

To quantify differences in protein accumulation due to heat 
and drought treatments, the barley leaf proteome of samples 





Fig. 2. Representative Coomassie-stained 12% SDS-PAGE containing 400 \ig total leaf protein. Numbered arrows indicate spots that 
were identified by MS and significantly regulated between conditions or between genotypes (Supplementary Table S4). 
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harvested three days after treatment start was visualized using 
DIGE, as seen in a representative image in Supplementary Fig. 
S2. Samples from this time point were chosen based on the 
RWC and PI measurements; 3 days after treatment start was 
the first time point where the combination treatment resulted 
in significant changes in RWC and PI in both genotypes com- 
pared to controls. Spot detection on a composite fluorescent 
image comprising all superimposed gel images resolved 1005 
distinct spots. Based on DIGE, this study identified 305 spots 
significantly differentially regulated by the heat treatment, 
473 spots different between genotypes, and 35 spots differ- 
ent due to the interaction between temperature and genotype. 
However, no spots were found to be significantly regulated by 
the drought treatment. Of the 305 spots found to be differen- 
tially regulated by heat, 99 were identified via mass spectrome- 
try (32 downregulated, 67 upregulated). Of the 473 spots found 
to be differentially regulated between genotypes, 125 were iden- 
tified (90 downregulated in Keel, 35 downregulated in Arta). 
Additionally, mass spectrometry helped identify 14 out of the 
35 proteins controlled by a significant interaction effect of the 
genotype and heat treatment (Supplementary Figs. S3 and S4). 

The biological functions of differentially regulated pro- 
teins included roles in metabolism, photosynthesis, transport, 
response to abiotic stimulus, and response to stress (Table 3 
and Supplementary Table S4). The proteomic analysis revealed 
the upregulation of structural components of the photosys- 
tem, i.e. the chlorophyll a-b binding protein of LHCII type III 
(Lhcb3) (spot 870), chloroplast oxygen-evolving enhancer pro- 
tein 1 (PsbO) (spots 846, 847, and 851) and oxygen-evolving 
enhancer protein 2 (PsbP) (spots 97 and 221). While the overall 
trend of genotype-specific protein regulation was to be down- 
regulation in Keel, photosynthesis-associated proteins had 
the opposite trend; 14 out of the 20 photosynthetic proteins 
were detected to be upregulated in Keel. In addition, proteins 
involved in detoxification, energy metabolism, and protein 
biosynthesis were differentially regulated by heat and between 
genotypes. These involved, for example, the upregulation of 
glycolytic proteins under heat, i.e. fructose-bisphosphate aldo- 
lase (spot 788) and cytosolic glyceraldehyde-3-phosphate dehy- 
drogenase (spots 327 and 330). In addition, chaperones (spots 
564 and 939), proteases (spots 586 and 527), elongation fac- 
tors (spots 533 and 723), and initiation factor 4A (spot 703 
and 710) were upregulated under heat. A large number of pro- 
teins differentially regulated between genotypes had functions 
in photorespiration and these were primarily upregulated in 
Arta as compared to Keel irrespective of the treatment. The 
proteomic data thus revealed a large number of proteins dif- 
ferentially regulated between genotypes and between heat and 
control conditions and the majority of these proteins had func- 
tions in photosynthesis, energy metabolism, and detoxification. 

Discussion 

Drought affected plant growth but not photosynthesis 

The reductions in grain yield due to the drought and heat 
treatments, while of similar magnitude, were due to different 



changes in yield component traits. Drought had major effects 
on plant growth, notably biomass and spike number, while 
the heat stress primarily affected the generative organs of the 
plant, the number of aborted spikes, and kernel weight. In 
addition, while drought did not significantly impact photo- 
synthetic efficiency, photosynthesis was significantly com- 
promised under heat and under the combination treatment. 
Consequently, drought affected plant growth but not pho- 
tosynthesis. The strong reduction in growth and simultane- 
ous maintenance of photosynthesis under drought as seen 
in the present study confirm recent studies in Arabidopsis, 
which demonstrated that mild drought primarily affected 
plant growth, but had minor effects on the photosynthesis 
rate (Muller et al, 2011; Skirycz et al, 2011; Verelst et al, 
2012). The results are thus consistent with the hypothesis that 
plants reduce their growth as a primary adaptation response 
to stress rather than as a secondary consequence of resource 
limitations (Muller et al, 201 1). In particular, under drought, 
biomass was positively correlated with yield. Thus, under 
non-lethal stress, limiting growth reduction might provide a 
strategy to increase productivity under stress. In the present 
study, drought also reduced the number of spikes presumably 
through reducing tiller number. The reduction in spike num- 
ber due to the inhibition of tillering is a known reaction to 
drought (El Soda et al, 2010), but its molecular basis is not 
established. However, drought-induced inhibition of tillering 
may share similar pathways to the inflorescence-induced inhi- 
bition of tillering. The repression of tillering that occurs dur- 
ing anthesis is thought to be due to a combination of auxin 
signalling and resource competition between the apical buds, 
which form new tillers, and the stem apex (Jewiss, 1972). 

Analysis of leaf proteomes under drought did not reveal 
significant changes of protein abundance as compared to 
control conditions. This is in contrast to the large number of 
transcripts differentially regulated in barley plants subjected 
to drought at the generative stage as seen by Guo et al. (2009). 
Homeostasis of the proteome could imply that the plants 
had acclimated to the drought stress (Harb et al, 2010). 
Acclimation to drought on the physiological level in the 
present study was evident in the maintenance of photosyn- 
thesis and water status under drought. Differences in stress 
symptoms might also be due to differences in the severity and 
duration of the stress application as compared to other stud- 
ies (Talame et al, 2007). However, the present data may also 
indicate a higher stability of the proteome as compared to 
the transcriptome under drought stress. Interestingly, a recent 
study has shown that while 1222 transcripts were differen- 
tially regulated in Arabidopsis subjected to mild drought, 34 
out of 208 1 proteins were significantly changed with an aver- 
age regulation factor of 1.5 (Baerenfaller et al., 2012). This 
suggests a higher stability of the proteome as compared to 
the transcriptome in response to environmental perturbation. 
These data suggest that barley adapted to non-lethal drought 
by avoidance mechanisms, in particular through the reduc- 
tion of growth. The resistance strategy ensured homeostasis 
of the cell which was reflected in the maintenance of photo- 
synthesis and stability of the proteome under drought. 
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Table 3. Proteins discussed in the text which were differentially regulated by temperature (T), genotype (G), or an interaction effect (GxT) 

Proteins quantified by DIGE, identified via mass spectrometry and grouped according to their biological function. Spot number (No.) in Fig. 2 
is given in addition to the Uniprot protein name and accession number. Predicted molecular weight (MW), isoelectric point (pi), Mascot score 
(Score), and percentage sequence coverage (SC%) are based on Mascot searches. The regulation factors, the log2 fold change in protein 
expression, are given for plants grown at 36 °C over plants grown at 21 °C (36/21) across both genotypes, for Keel plants over Arta plants (K/A) 
across all treatments, for heat-treated Keel plants over heat-treated Arta plants (K36/A36) and for control Keel plants over control Arta plants 
(K21/A21). Regulation factors corresponding to significant (P < 0.05) changes in expression are underlined. The complete list of differentially 
regulated proteins can be found in Supplementary Table S3. 



No. ■ 



Protein name 



UniReflOO MW 



pi Score SC% SE 36/21 



K/A 



K36/A36 K21/A21 



Carbohydrate metabolic process 



327 



330 



788 



Glyceraldehyde-3-phosphate 

dehydrogenase, 

cytosolic 

Glyceraldehyde-3-phosphate 

dehydrogenase, 

cytosolic 

Fructose- 

bisphosphatealdolase 



Photosynthesis 

97 Oxygen-evolving 

enhancer protein 2, chloroplastic 
221 Oxygen-evolving 

enhancer protein 2, chloroplastic 

846 Predicted protein 

847 Chloroplast oxygen-evolving 
enhancer protein 1 

851 Predicted protein 

870 Chlorophyll a-b binding protein 

of LHCII type III, chloroplastic 
Response to abiotic stimulus 



586 



527 



564 
939 

Translation 

533 
703 

710 
723 



ATP-dependent zinc metallopro- 
tease FTSH 1 , chloroplastic 
ATP-dependent Clp protease 
ATP-binding subunit clpA homo- 
logue CD4B, chloroplastic 
Chloroplast heat shock 
protein 70 
Heat-shock protein 

Elongation factor EF-G 
Eukaryotic initiation 
factor 4A 

Eukaryotic initiation factor 4A 
Elongation factor Tu 



P26517 

P26517 

F2ELD1 

Q00434 

Q00434 

F2CRK1 
A5JV93 

F2CRK1 
P27523 

Q5Z974 
P31542 

A4ZYQ0 

Q43638 

Q9SI75 
P41378 

P41378 
Q8W2C3 



122.3 7.0 154.6 2.4 



48.3 7.0 49.0 2.6 



41.9 7.5 304.6 10.! 



96.0 

96.0 

34.4 
68.7 

34.4 
91.6 

80.0 
173.3 

134.0 

103.9 

61.3 
143.9 

68.3 
144.5 



6.9 
6.9 



6.9 



44.9 0.9 
91.5 7.0 



5.6 306.6 11.6 
7.0 198.6 3.7 



5.6 199.0 11.6 



7.0 108.6 



6.9 121.0 



1.2 



6.9 74.2 9.2 

6.9 124.0 7.9 

7.0 170.2 1.8 

7.0 122.8 2.2 



64.2 16.0 



7.0 



74.4 



6.9 298.2 



4.8 



2.0 
2.4 



1.29 



T/G/GxT 1. 



T/G 

T/G 

T/G 
T/G 

T/G 
T/G 



T/G 
T 

T/G 



1.22 

1.46 

1.29 

1.29 
1.55 

1.63 
1.75 

2.69 
-2.07 

1.85 
2.90 



T/G/GxT 2.41 
T/G/GxT 1.73 



1.02 

-1.97 

1.09 

2.14 

1.86 

1.80 
2.53 

2.39 
1.98 



G 
G 



-1.11 
-1 .08 



-1.10 

-1 .43 

-1 .45 
-1.65 

-1 .56 
-1 .34 



-1.18 

-2.36 

1.07 

2.10 

1.84 

1.82 
2.52 

2.28 
2.05 



-1 .29 -1 .37 
-1.76 -2.35 



-1.21 

-1 .54 

-1.52 
-1.96 

-1 .86 
-1 .40 



1.30 



-1.49 



1.10 



2.21 



1.78 
2.55 

2.58 
1.86 



-1.12 
-1.55 

1.08 
-1.16 



-1.29 
-1.24 



-1 .34 
-1.28 



Heat affected photosynthesis rates 

The heat treatment caused significant changes in reproductive 
organs and affected TKW and spike abortion rates, while the 
vegetative traits biomass and plant height remained relatively 
unaffected. At the same time, heat caused a significant reduc- 
tion of RWC and photosynthetic efficiency. The RWC was likely 
reduced by an increase in transpiration as seen in the reduced 
leaf temperature under heat as compared to the ambient tem- 
perature. As ambient temperature increases, plants attempt to 
cool themselves by opening stomata and increasing transpira- 
tion (Schulze et al., 1973). However, the RWC under heat was 
not significantly different from the RWC under drought con- 
ditions, while photosynthesis rates were significantly different 



between these two stresses. This suggested that photosynthesis 
was significantly reduced under heat not because of suboptimal 
water content in the cells but rather by a direct detrimental effect 
of heat on the photosynthetic apparatus. These findings are in 
line with previous observations that the optimum temperature 
for photosynthesis in barley is 20 °C and photosynthesis rates 
are decreased by more than 50% at 35 °C (Todd, 1982). The det- 
rimental effect of heat on the photosynthetic apparatus was also 
indicated by the proteome data, which demonstrated an increase 
in the turnover of photosynthesis-related proteins under heat. 

A reduction in Fv/Fm and PI in the leaves and senescence of 
the lower leaves (Fig. 1) indicated that the ability to photosyn- 
thesize was permanently decreased under heat. In addition, 
the PI was significantly reduced under heat which suggests that 
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both the light-dependent and light-independent mechanisms, 
e.g. carbon fixation, were damaged or inhibited due to high 
temperature. Thus, the heat treatment had a detectable effect 
on photoinhibition and the availability of carbon dioxide was 
probably limited by stomata closure. Evidence of a reduction 
of the carbohydrate pool in heat-treated plants was also pro- 
vided by the reduction in kernel weight and by the upregula- 
tion of glycolysis enzymes observed in the proteomic data. 

Heat also reduced total grain yield by increasing the abor- 
tion rate of spikes. Floret fertility and grain setting, as meas- 
ured by the number of aborted spikes, was reduced under the 
heat treatment. Reproductive growth is known to be more 
sensitive to heat stress than vegetative growth in barley. In 
particular, anthers are prone to growth inhibition and lose 
the ability to produce pollen when heat stressed (Oshino 
et al., 2007). Heat thus caused yield reductions primarily by 
affecting photosynthesis, carbon fixation, and floret fertility. 
Strong effects of the short-term heat treatment on physiology 
indicated that the homeostasis of the cells was disturbed. 

As in stress-prone environments drought and heat often co- 
occur, this study compared the combined effects of drought 
and heat with that of the single stresses on plant performance. 
The changes in plant performance were greater under the 
combination treatment than under heat or drought alone 
(Table 2). The strongest effect of the combination treatment 
was observed for PI suggesting that photosynthesis was most 
sensitive to the simultaneous application of drought and heat, 
presumably because leaf cooling through transpiration was 
severely reduced by the limited water supply. Reduced tran- 
spiration and leaf cooling in the combination treatment was 
also suggested by the difference in leaf temperature between 
the heat (33-34 °C) and the combination treatment (36-37 °C) 
(Table 2). In addition, grain yield was significantly more 
reduced by the combination treatment than by the single 
stresses, presumably because drought and heat affected unique 
yield component traits such as spike number and spike fertility, 
which resulted in an overall reduced yield in the combination 
treatment. 

Proteomic basis of morphological plasticity and 
physiological responses to heat stress 

The proteome analysis identified significant differences in 
protein abundance only under heat, and not under drought. 
This reflects the strong effects of heat on plant physiology in 
contrast to the physiological homeostasis seen under drought. 
Under heat, a large number, 99 out of 296 detected proteins 
were differentially regulated, and these had predominantly 
functions in photosynthesis, detoxification, energy metabo- 
lism, and protein biosynthesis (Table 3 and Supplementary 
Table S4). In many cases, several spots were detected per pro- 
tein. These may represent close homologues, which could not 
be resolved based on the mass spectrometric data, but may 
also be isoforms due to differential post-translational modi- 
fications (Rohrig et al., 2006). The high frequency of differ- 
entially regulated proteins with functions in photosynthesis 
may be explained by the enrichment of plastid proteins in the 
present leaf proteome analysis (Supplementary Fig. S2). The 



proteomic analysis revealed structural components of the 
light-harvesting complex (Lhcb3) and the oxygen-evolving 
complexes (PsbO and PsbP) as being significantly upregu- 
lated under the heat treatment. It is tempting to speculate 
that the upregulation of Lhcb3, PsbO, and PsbP might be 
due to the de novo synthesis of peptides that are en route to 
replace them and the lag in degradation of damaged proteins 
removed from the photosystem. The upregulation of Lhcb3, 
PsbO, and PsbP proteins was therefore likely due to their 
increased replacement after damage by heat as indicated by 
the chlorophyll fluorescence data. The reduction of the PI 
due to the heat treatment was indicative of inhibition of the 
light-independent reactions of photosynthesis. Rubisco, the 
enzyme catalysing the rate-limiting step of carbon fixation, is 
inhibited by side-products of photorespiration which stabilize 
the active site of Rubisco in a closed conformation (Pearce 
and Andrews, 2003). Interestingly, Rubisco activase B showed 
the highest upregulation under heat, while Rubisco activase 
A was downregulated under heat (Table 3). Rubisco activase 
frees the catalytic site of Rubisco from inhibitory sugar phos- 
phates by forcing the active site into an open conformation 
(Portis et al, 2007). Rubisco activases have been characterized 
on the genomic level in barley (Rundle and Zielinski, 1991), 
wheat (Law and Crafts-Brandner, 2001), rice (Orysa sativa L.) 
(To et al, 1999), cotton (Gossypium hirsutum L.) (Salvucci 
et al., 2003), maize (Zea mays L.) (Ayala-Ochoa et al., 2004), 
and Arabidopsis (Werneke et al, 1989). Rubisco activase 
A with two different splice variants is present in all mentioned 
species and is known to be heat inactivated starting at 35 °C 
(Crafts-Brandner et al, 1997). In contrast, Rubisco activase 
B has only been detected in barley, wheat, maize, and cotton 
but not in rice or Arabidopsis. It has already been shown that 
Rubisco activase B is induced by heat on the transcript level 
in wheat (Wang et al, 201 1) and on the protein level in cotton 
(Law et al., 2001). However, the thermostability of Rubisco 
activase B is currently untested. The upregulation of Rubisco 
activase B under heat treatment suggested a specific role for 
Rubisco activase B in maintaining the activity of Rubisco 
under high temperature conditions, possibly by being more 
thermostable than Rubisco activase A. 

Interestingly, the proteomic data showed that several gly- 
colytic proteins were upregulated under the heat treatment. 
The prime functions of glycolysis are to generate carbon skel- 
etons, reductants, and ATP, which can confer a bioenergetic 
advantage that can extend the survival time of plant cells that 
have become ATP-depleted due to environmental stresses. 
Increased levels of GAPDH transcripts have been observed 
under environmental stress conditions, such as dehydration 
in Cratero stigma plantagineum (Velasco et al., 1994), during 
heat shock in Arabidopsis plants (Yang et al., 1993), and dur- 
ing anaerobic stress in maize (Chang et al., 2000). Thus, there 
is supporting evidence for the activation of ATP-generating 
pathways under different stresses, presumably to cope with a 
higher demand for ATP to maintain homeostasis under stress 
conditions. A higher energy demand was also suggested by 
the upregulation of ATP synthase subunit alpha in the mito- 
chondria and plastids upon heat stress (Table 3). For example, 
higher transcript levels of ATP synthase were observed in rice 
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upon salt and osmotic stress (Zhang et al., 2006). A higher 
energy was probably required for a higher rate of protein deg- 
radation and biosynthesis as suggested by the upregulation 
of chaperones, proteases, elongation factors, and initiation 
factor 4A. Analysis of the proteome of heat-stressed barley 
plants thus suggested thermolability and a more rapid protein 
turnover of photosynthesis-related proteins, while drought 
did not affect cellular homeostasis. 

Barley genotypes Arta and Keel responded uniquely to 
heat and drought treatments 

The identification of genetic differences in stress responses 
is an important basis for improving plant performance 
under stress. The differences in responses to drought and 
heat, in particular in growth, spike fertility, and chloro- 
phyll fluorescence seen between Arta and Keel, suggested 
that the two genotypes have unique mechanisms for cop- 
ing with environmental stresses. Arta was characterized 
by a stronger growth reduction under drought than Keel. 
Under heat, Arta showed a higher spike abortion rate as 
compared to Keel. Finally, under the combination treat- 
ment Arta had a higher Fv/Fm than Keel. Genetic variation 
in the efficiency of photosynthesis under stress has already 
been detected in diverse barley genotypes (Oukarroum 
et al., 2007; Ahmed et al., 2013). Interestingly, Keel tended 
to show higher yield under drought and heat alone, while 
Arta tended to show a higher yield under the combination 
treatment. Differences in growth and yield under the dif- 
ferent stress regimes suggested that Keel is better adapted 
to mild stress as it maintains growth and yield. In contrast, 
under more severe stress as represented by the combina- 
tion treatment Arta performed better presumably through 
a stronger morphological adaptation and maintenance of 
photosynthesis. The genetic difference was supported by the 
high number of proteins differentially regulated between 
the genotypes. Interestingly, a large number of these dif- 
ferential proteins had roles in photorespiration, and these 
were primarily downregulated in Keel. Especially under 
stress conditions that lead to reduced rates of photosyn- 
thetic C0 2 assimilation, photorespiration serves as energy 
sink preventing the over-reduction of the photosynthetic 
electron transport chain and photoinhibition. Differences 
in the efficiency of photorespiration between Arta and Keel 
may explain the difference in photosynthetic performance 
under the combination stress. This study also observed 14 
protein spots regulated by interaction effects between geno- 
type and environment which are thus potential outputs of 
unique adaptations to heat stress that have evolved between 
Arta and Keel (Supplementary Figs. S3 and S4). Arta and 
Keel are parents of a recombinant inbred population, which 
shows transgressive segregation for agronomic performance 
under drought in Mediterranean environments (J. Rollins, 
B. Drosse M. A. Mulki, S. Grando, M. Baum, M. Singh, S. 
Ceccarelli, M. von Korff, in revision) unpublished results). 
Identification of the unique stress adaptation strategies in 
Arta and Keel will allow further dissection of the genetic 
basis of this transgressive performance in the offspring. 



Conclusion 

The data presented in this study suggested that barley has 
adapted to non-lethal drought by avoidance mechanisms, such 
as the reduction of growth which allowed the plants to main- 
tain a cellular homeostasis as seen in the stability of photosyn- 
thesis and of the proteome under drought. In contrast, heat 
affected RWC, photosynthesis, and floret fertility and caused 
a rapid turnover of photosynthesis-related proteins. Based on 
the protein changes observed, it is proposed that important 
heat-tolerance mechanisms include protein quality control 
and de novo synthesis, which cause a higher energy demand as 
seen in the upregulation of ATP generating pathways. Reduced 
C0 2 availability due to stomatal closure under heat was coun- 
ter-balanced by the activation of Rubisco by Rubisco activase 
B. Inhibition of the light-independent reactions of photosyn- 
thesis likely caused the production of radical oxygen species 
as seen in the increase of scavenging and detoxifying enzymes. 
Genetic variation in stress responses between Arta and Keel, in 
particular in growth, spike fertility, and photosynthesis can be 
exploited in future crop breeding efforts. 
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